Recombination of Shower Partons in Fragmentation Processes 
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We develop the notion of shower partons and determine their distributions in jets in the framework 
of the recombination model. The shower parton distributions obtained render a good fit of the 
fragmentation functions. We then illustrate the usefulness of the distributions in a problem where a 
jet is produced in the environment of thermal partons as in heavy-ion collisions. The recombination 
of shower and thermal partons is shown to be more important than fragmentation. Application to 
the study of two-particle correlation in a jet is also carried out. 
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I. INTRODUCTION 

The theoretical description of hadron production at 
large transverse momentum (pr) in either hadronic or 
nuclear collisions at high energies is traditionally framed 
in a two-step process that involves first a hard scattering 
of partons, followed by the fragmentation of the scat- 
tered parton to the detected hadron 0, 0|. The first 
part is calculable in perturbation QCD, while the second 
part makes use of fragmentation functions that are de- 
termined phenomenologically. Such a production mech- 
anism has recently been found to be inadequate for the 
production ofparticles at intermediate pt in heavy-ion 
collisions 0. 0- Q- Instead of fragmentation it is the 
recombination of partons that is shown to be the more 
appropriate hadronization process, especially when the 
soft partons are involved. Although at extremely high 
Pt fragmentation is still dominant, it is desirable to have 
a universal description that can be applied to any pt, 
based on the same hadronization scheme. To achieve 
that goal it is necessary that the fragmentation process 
can be treated as the result of recombination of shower 
partons in a jet. The purpose of this paper is to take 
that first step, namely: to introduce the notion of shower 
partons and to determine their distributions in order to 
represent the phenomcnological fragmentation functions 
in terms of recombination. 

The subject matter of this work is primarily of interest 
only to high-energy nuclear collisions because hadroniza- 
tion in such processes is always in the environment of soft 
partons. Semi- hard shower partons initiated by a hard 
parton can either recombine among themselves or recom- 
bine with soft partons in the environment. In the former 
case the fragmentation function is reproduced, and noth- 
ing new is achieved. It is in the latter case that a very new 
component emerges in heavy-ion collisions, one that has 
escaped theoretical attention thus far. It should be an 
important hadronization process in the intermediate pr 
region. Our main objective here is to quantify the prop- 
erties of shower partons and to illustrate the importance 
of their recombination with thermal partons. The actual 
application of the shower parton distributions (SPD) de- 



veloped here to heavy-ion collisions will be considered 
elsewhere 

The concept of shower partrons is not new, since at- 
tempts have been made to generate such partons in 
pQCD processes as far as is permitted by the validity of 
the procedure. Two notable examples of such attempts 
are the work of Marchesini and Webber and Geiger 
0. However, since pQCD cannot be used down to the 
hadronization scale, the branching or cascading processes 
terminate at the formation of color-singlet pre-hadronic 
clusters, which cannot easily be related to our shower par- 
tons and their hadronization. We shall discuss in more 
detail at the end of Sees. Ill and IV the similarities and 
differences in the various approaches. 



II. RECOMBINATION MODEL FOR 
FRAGMENTATION 

The fragmentation of a parton to a hadron is not a 
process that can be calculated in pQCD, although the 
Q 2 evolution of the fragmentation function (FF) is calcu- 
lable. The FF's are usually parameterized by fitting the 
data from e + e~ annihilations 0, 0, ^| as well as from 
pp and e ± p collisions Although the QCD processes 
of generating a parton shower by gluon radiation and 
pair creation cannot be tracked by perturbative meth- 
ods down to low virtuality, we can determine the SPD's 
phenomenologically in much the same way that the FF's 
themselves are, except that we fit the FF's, whereas the 
FF's are determined by fitting the data. An important 
difference is that both the shower partons and their dis- 
tributions are defined in the context of the recombination 
model, which is the key link between the shower partons 
(inside the black box called FF) and the observed hadron 
(outside the black box). 

In the recombination model the generic formula for a 
hadronization process is 01 

f x dx\ f x dx2 

xD[X) = / / F q qi{Xi,X2)R(Xi,X2,X) , (1) 

Jo x l Jo x 2 

where F qq -i(xi, X2) is the joint distribution of a quark q 
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at momentum fraction x\ and an antiquark q' at x 2 , and 
R(xi, x 2 , x) is the recombination function (RF) for the 
formation of a meson at x. We have written the LHS of 
Eq. Q as xD(x), the invariant FF, but the RHS would 
have the same form if the equation were written for the 
inclusive distribution, xdN/dx, of a meson produced in a 
collisional process. In the former case of fragmentation, 
F q qi refers to the shower partons initiated by a hard par- 
ton. In the latter case of inclusive production, F qq i refers 
to the q and q' that are produced by the collision and are 
to recombine in forming the detected meson. The equa- 
tions for the two cases are similar because the physics of 
recombination is the same. In either case the major task 
is in the determination of the distribution F qq i . 

We now focus on the fragmentation problem and re- 
gard Eq. Q as the basis of the recombination model for 
fragmentation. The LHS is the FF, known from the pa- 
rameterization that fits the data. The RHS has the RF 
that is known from previous studies of the recombination 
model 0, 0| and will be specified in the next section. 
Thus it is possible to determine the properties of F q qi 
from Eq. To facilitate that determination we shall 
assume that F qq i is factorizable except for kinematic con- 
straints, i.e., in schematic form we write it as 

F$(x 1 ,x 2 ) = Sf(x 1 )S?(x 3 ) , (2) 

where Sf{x\) denotes the distribution of shower parton 
q with momentum fraction xi in a shower initiated by 
a hard parton i. The exact form with proper kinematic 
constraints will be described in detail in the next section. 
Here we remark on the general implications of Eqs. 
and {2J. 

The important point to emphasize is that we are intro- 
ducing the notion of shower partons and their momentum 
distributions Sf(xi). The significance of the SPD is not 
to be found in problems that involve only the collisions 
of leptons and hadrons, for which the fragmentation of 
partons is known to be an adequate approach, and the re- 
combination of shower partons merely reproduces what 
is already known. The knowledge about the SPD be- 
comes crucial when the shower partons recombine with 
other partons that are not in the jet but are in the ambi- 
ent environment. We shall illustrate this important point 
later. 

It should be recognized that the SPD that we shall de- 
termine through the use of Eqs. {TJ and {2Jl depends on 
the specific form of R(xi, x 2 , x), which in turn depends 
on the wave function of the meson produced. It would 
be inconsistent to use our Sf given below in conjunction 
with some approximation of the RF that differs signif- 
icantly from our R. The recombination of two shower 
partons must recover the FF from which the SPD's are 
obtained. 

Finally, we remark that Sf should in principle depend 
on Q 2 at which the D(x,Q 2 ) is used for its determina- 
tion, since Q 2 evolution affects both. It is outside the 
scope of this paper to treat the Q 2 dependence of Sf. 



Our aim here is to show how Sf can be determined phe- 
nomenologically, and how it can be applied, when Q 2 is 
fixed. The same method can be used to determine Sf at 
other values of Q 2 . In practice, the Q 2 dependence of 
Sf is not as important as the inclusion of the role of the 
shower partons in the first place at any reasonably ap- 
proximate Q 2 in heavy-ion collisions where hard partons 
are produced in a range of transverse momentum. 

III. SHOWER PARTON DISTRIBUTIONS 

In order to solve Eqs. {TJ and {2J for Sf, we first point 
out that there are various D(x) functions corresponding 
to various fragmentation processes. We shall select five 
of them, from which we can determine five SPD's. Three 
of them form a closed set that involves no strange quarks 
or mesons. Let us start with those three. Consider the 
light quarks u, d, u, d, and gluon g. They can all frag- 
ment into pions. To reduce them to three essential FF's, 
we consider the three basic types D v , Dg and Dq, that 
correspond to valence, sea and gluon fragmentation, re- 
spectively. If the fragmenting quark has the same flavor 
as that of a valence quark in n, then the valence part 
of the fragmentation is described by Dy, e.g., u v — > 7r + , 
d v — * 7r~, u v — > 7r~ , the subscript v referring to the va- 
lence component. All other cases of quark fragmentation 
are described by Dg, e.g., u — > 7r~, d — ► 7r + , u — ► 7r + . 
If the initiating parton is a gluon, then we have Dq for 
any state of ir. Those FF's are given by Ref. 9] in para- 
metric form. We shall use them even th oug h they are 
older than the more recent ones |ToL ITTl Il4j . which do 
not give the D v and Dg explicitly. Our emphasis here is 
not on accuracy, but on the feasibility of extracting the 
SPD's from the FF's of the type discussed above. We 
shall determine Sf from the BKK parameterization 
with Q 2 fixed at 100 GeV 2 and demonstrate that the use 
of shower partons is important in heavy-ion collisions. 

For the SPD's we use the notation K^s and L for 
valence and sea-quark distributions, respectively, in a 
shower initiated by a quark or antiquark, and G for any 
light quark distribution in a gluon-initiated shower. That 
is, for example, K NS = S™ v , L = S*, G = S q . It should 
be recognized that L also describes the sea quarks of the 
same flavor, such as S™""" , so that the overall distribution 
of shower quark that has the same flavor as the initiating 
quark (e.g. u — > u) is given by 

K = K NS + L . (3) 

It is evident from the above discussion that there is 
a closed relationship that is independent of other un- 
knowns. It follows from Eq. when restricted to sea- 
quark fragmentation: 

xD s {x) = f^l^lL(x 1 )L(-^-)R 7r (x u x 2 ,x) .(4) 
J Xi x 2 \l-x 1 J 

The sea-SPD L(z) can be determined from this equation 
alone. In Eq. {3J we have exhibited the argument of the 
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second L function that reflects the momentum constraint, 
i.e., if one shower parton has momentum fraction xi, then 
the momentum fraction of the other recombining shower 
parton cannot exceed 1 — x\, and can only be a fraction 
of the balance £2/(1 — x\). Symmetrization of x\ and x 2 
is automatic by virtue of the invariance of R 7T (x± 1 X2, x) 
under the exchange of x\ and X2 ■ 

After L(z) is determined from Eq. Q, we next can 

I 



{K NS (x 1 ),L(x 2 )} = - 



Finally, we have the closed equation for the gluon- 
initiated shower 

xD G {x) = { ^1^1 G (x 1 )g(-^-)r 7T (x 1 ,x 2 ,x).(7) 
J xi x 2 \l-xij 

In this non-strange sector we have 3 SPD's (L, K^s and 
G) to be determined from the 3 phenomenological FF's 
(Dg, D v and D G ). 

In extending the consideration to the strange sector, 
we must make use of L and G determined in the above 
set and two new FF's Df and Dq that describe the frag- 
mentation of a non-strange quark and gluon, respectively, 
to a kaon. That is, we have 

xDg(x) = [ ^Hl^Hl {L(xi),L s (x 2 )} Rk(xi,x 2 ,x), (8) 
J xi x 2 

xD«(x)= [ ^^-{G(x 1 ),G s (x 2 )}R K (x 1 ,x 2 ,x),(9) 
J xi x 2 

where L s and G s are two additional SPD's specifying the 
strange quark distributions in showers initiated by non- 
strange and gluon partons, respectively. Rk is the RF 
for kaon. 

To complete the description of the integral equations, 
we now specify the RF's. They depend on the square 
of the wave functions of the mesons, 7r and K, whose 
structures in momentum space have been quantified in 
the valon model [l2L . Unlike the case of the proton, 
whose structure is well studied by deep inelastic scatter- 
ing so that the valon distribution can be obtained from 
the parton distribution functions 0] , the RF for the pion 
relies on the parton distribution of the pion probed by 
Drell-Yan process [l^. The derivation of the RF's for 
both 7r and K is given in [l^; they are 

R K (xi,x 2 ,x) = — —51 1 1 , (10) 

X z \ X X I 



obtain K^s from 

f dx dx 

xD v (x) = / — - — - {K NS (xi),L(x 2 )} R 7r (x 1 ,x 2 ,x),(5) 

J Xi x 2 

where the curly brackets define the symmetrization of the 
leading parton momentum 



(6) 

I 

The 5 functions guarantee the momentum conservation 
of the recombining quarks and antiquarks, which are 
dressed and become the valons of the produced hadrons. 

Since the recombination process involves the quarks 
and antiquarks, one may question the fate of the gluons. 
This problem has been treated in the formulation of the 
recombination model [T^ , where gluons are converted to 
quark-antiquark pairs in the sea before hadronization. 
That is, the sea is saturated by the conversion to carry 
all the momentum, save the valence parton. Such a pro- 
cedure has been shown to give the correct normalization 
of the inclusive cross section of hadronic collisions [l2| . 
In the present problem of parton fragmentation we imple- 
ment the recombination process in the same framework, 
although gluon conversion is done only implicitly. What 
is explicit is that the gluon degree of freedom is not in- 
cluded in the list of shower partons. It means that in 
the equations for D v , Dg, and D G (and likewise in the 
strange sector) only K^s, L and G appear; they are the 
SPD's of quarks and antiquarks that are to recombine. 
Those quarks and antiquarks must include the converted 
sea, since they are responsible for reproducing the FF's 
through Eqs. 10} , (J5J and J7J without gluons. Thus the 
shower partons whose momentum distributions we cal- 
culate are defined by those equations that have no gluon 
component for recombination, and would not be the same 
as what one would conceptually get (if it were possible) 
in a pQCD calculation that inevitably has both quarks 
and gluons. 

It should be noted that our procedure of converting 
gluons to qq pairs is essentially the same as what is done 
in p| , whose branching processes terminate at the thresh- 
old of the non-perturbative regime. In that approach 
nearby quarks and antiquarks that are the products of 
the conversion from different gluons form color-singlet 
clusters of various invariant masses that subsequently de- 
cay (or fragment as in strings) sequentially through res- 
onances to the lowest lying hadron states ^tJ- Similar 
but not identical approach is taken in [l8| , where gluons 
are not directly converted to qq pairs, but are either ab- 
sorbed or annihilated byg + g—>q + q Born-diagram 



K N s{x\)L 



x 2 

1 — X\ 



Kns 



X\ 

1 - x 2 



L(x 2 ) 



4 



processes. 

IV. RESULTS 

We now proceed to solve the integral equations for the 
five FF's, which are known from Ref. [!|. Those equa- 
tions relate them to the five unknown SPD's: K^s, L, 
G, L s and G s . If those equations were algebraic, we obvi- 
ously could solve them for the unknowns. Being integral 
equations, they can nevertheless be "solved" by a fitting 
procedure that should not be regarded as being unsat- 
isfactory for lack of mathematical rigor, since the FF's 
themselves are obtained by fitting the experimental data 
in some similar manner. Indeed, the FF's in the next-to- 
leading order are given in parameterized forms Q 

D%(x) = Nx a (l-xf(l + xy (12) 

where the parameters for Q 2 — 100 GeV 2 are given in 
Table I for k = S, V, G and h = it, K. 

TABLE I: Parameters in Eq. O for Q 2 = lOOGeV 2 . 





N 


a 




7 


Ds 


2.7236 


-0.734 


3.384 


-5.471 


Dl 


0.2898 


-1.040 


1.608 


-0.111 


Dh 


0.7345 


-1.112 


2.547 


-0.541 


D§ 


0.2106 


-1.005 


2.548 


-0.620 


D% 


0.0768 


-1.481 


2.489 


-0.778 
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FIG. 1: Fragmentation functions, as parametrized in are 
shown in symbols, while those calculated in the recombination 
model are shown by the solid lines. All curves are for Q 2 — 
lOOGeV 2 . 

All five SPD's are denoted collectively by Sf with i = 
q, q, g and j = q, s, q, s, where q can be either u or d. If 
in i the initiating hard parton is an s quark, it is treated 



as q. That is not the case if s is in the produced shower. 
Our parametrization of Sf has the form 

Sl{z) = Az a {l-z) b {l + cz d ), (13) 

where the dependences of the parameters A, a, etc. on 
i and j are not exhibited explicitly, just as in Eq. I|12|l . 
Substituting Eqs. ljTU|) and (fTTft into |gJ|-@, we can de- 
termine the parameters one equation at a time, i.e., L 
from Dg, G from D^, and then K^s from Dy and so 
on. In most cases it can be shown that the x — > 1 limit 
requires b — (3. The final results of the fits are shown in 
Fig. 1 with the corresponding parameters given in Table 
II. 



TABLE II: Parameters in Eq. (13 . 





A 


a 


b 


c 


d 


Kns 


0.333 


0.45 


2.1 


5.0 


0.5 


L 


1.881 


0.133 


3.384 


-0.991 


0.31 


G 


0.811 


-0.056 


2.547 


-0.176 


1.2 


L s 


0.118 


-0.138 


2.3 


0.90 


0.1 


Gs 


0.069 


-0.425 


2.489 


-0.5 


1.1 




FIG. 2: Shower parton distributions determined in the re- 
combination model, corresponding to the parameterization 
given in Table II. 

It is evident from Fig. 1 that all the fits are very good, 
except in the low x region of D v {x). In the latter case 
we are constrained by the condition 

f dz 

J ^fK NS {z) = 1 (14) 

that is imposed by the requirement that there can be 
only one valence quark in the shower partons. However, 
the fit for x > 0.4 is excellent, and that is the important 
region for the determination of Kns(z)- In application to 
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u — > 7r + , say, the u quark in the shower must have both 
valence and sea quarks so the shower distribution for the 
u quark is always the sum: K (z) = Kms{ z ) +L(z). Since 
L(z) is large at small z, and is accurately determined, the 
net result for K(z) should be quite satisfactory. 

It is remarkable how well the FF's in Fig. 1 are repro- 
duced in the recombination model. The corresponding 
SPD's that make possible the good fit are shown in Fig. 
2. They have very reasonable properties, namely: (a) va- 
lence quark is harder (b) sea quarks are softer, (c) gluon 
jet has higher density of shower partons, and (d) the den- 
sity of produced s quarks is lower than that of the light 
quarks. 
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FIG. 3: Parton distributions in transverse momentum kr for 
valence+sea quark (solid line), sea quark (dash-dot line) and 
thermal partons (dashed line). 

It is appropriate at this point to relate our approach 
to those of Marchesini- Webber Q and Geiger [8( , which 
are serious attempts to incorporate the QCD dynamics in 
their description of the branching and collision processes. 
The former is done in the momentum space only, whereas 
the latter is formulated in space-time as well as in mo- 
mentum space. The parton cascade model of Geiger is a 
very ambitious program that treats a large variety of pro- 
cesses ranging from e + e~ annihilation 18] to deep inelas- 
tic scattering 0] to hadronic and nuclear collisions [2(J . 
The evolution of partons is tracked by use of relativistic 
transport equations with gain and loss terms. Cluster 
formation takes into account the invariant distance be- 
tween near-neighbor partons. Cluster decay makes use 
of the Hagedorn spectrum and the particle data table. 
Because of the complexity of the problems both QCD 
models are implemented by Monte Carlo codes. The pre- 
dictive power of the models is exhibited as numerical out- 
puts that cannot easily be adapted for comparison with 
our results on the SPD's. Our approach makes no at- 
tempt to treat the QCD dynamics; however, the SPD's 
obtained are guaranteed to reproduce the FF's on the 
one hand, and are conveniently parameterized for use in 
other context that goes beyond fragmentation, as we shall 



show in the next section. From the way the color-singlet 
clusters are treated in the QCD models, it is clear that 
our shower partons do not correspond to the partons of 
those models at the end of their evolution processes, ex- 
cept in the special case when the cluster consists of only 
one particle. In our approach the non-perturbative part 
of how the shower partons dress themselves and recom- 
bine to form hadrons with the proper momentum-fraction 
distributions is contained in the RF's. Such shower par- 
tons that are ready to hadronize are sufficiently far from 
other shower partons as to be independent from them. 
In general, they cannot be identified with the q and q 
that form the color-singlet clusters in the QCD models, 
but are more closely related to the constituents of the fi- 
nal hadrons, as in the case of quarkonium formation |2l| . 
The distribution of those constituents in a hard-parton 
shower cannot be displayed in the QCD models, but are 
determined by us by solving Eqs. (@J and 10 



V. APPLICATIONS 

As we have stated in the introduction, the purpose 
of determining the SPD's is for their application to 
problems where the FF's are insufficient to describe the 
physics involved. We consider in this section two such 
problems as illustrations of the usefulness of the SPD's. 
The first is when a hard parton is produced in the en- 
vironment of thermal partons, as in heavy-ion collisions. 
The second is the determination of two-pion distribution 
in a jet. 

Let us suppose that a u quark is produced at kx = 10 
GeV/c in a background of thermal partons whose invari- 
ant kx distribution is 



T(kx) — kr— — = Cfc^e" 
akr 



k T /T 



(15) 



Let the parameters C and T be chosen to correspond to 
a typical situation in Au+Au collisions at y/s = 200 GeV 

a 



C = 23.2 GeV" 



T = 0.317 GeV. 



(16) 



The high-^T u quark generates a shower of partons 
with various flavors. Consider specifically u and d in 
that shower. The valence quark distribution is given by 
Kns{xi), while the d sea-quark distribution (including 
the ones converted from the gluons) is given by L(xi). 
In Fig. 3 we plot dN/krdkT for (a) u quark (valence and 
sea) in solid line, (b) d sea antiquark in dash-dot line, and 
(c) d thermal antiquark in dashed line. They correspond 
to k^ 2 x (invariant distributions K — Kns + L, L, and 
T, respectively), in which K(xi) and L{x\) are evalu- 
ated at k T = x x k^, with fc5? ax = 10 GeV/c. Note that 
the thermal distribution is higher than the shower parton 
distributions for kx < 1 GeV/c. That makes a crucial 
difference in the recombination of those partons. Such a 
thermal distribution is absent in pp collisions, whose soft 
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partons are at least two orders of magnitude lower. In 
e + e~ annihilation there are, of course, no soft partons at 
all. 



K and L from the D function in the first place, Eq. i|18|) 
can more directly be identified with 
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FIG. 4: Distributions of ir + in pr arising from thermal- 
shower recombination (solid line) and shower-shower recom- 
bination, i.e. fragmentation (dash-dot line). 

We now calculate the production of ir + from the as- 
semblage of u and d partons. The thermal-shower (TS) 
recombination gives rise to 

= — / dk T K(k T /k% ax )T( PT - br), (17) 
p T dp T P T Jo 

where Eq. 1)11) [1 has been used in an equation such as Eq. 
||TJ for xdN-x/dx, but expressed for dN^/pxdpT- Using 
Eqs. ©, ijBll and the parametrizations given in Table 
II, the integral in Eq. I|17|) can readily be evaluated. The 
result is shown by the solid line in Fig. 4. It is to be com- 
pared with the pt distribution from the fragmentation of 
the u quark to 7r + , which is 



dN s f 

7TJ_ 

PTdpT 



Pt -Jo 



dk T <K 



^,max 



p T 



a«) 



p T dpj 



/ 7 max\— 1 

{p T k T ) 



Pt 

Lmax 

tvrp 



Dl I 

S \ (.max 



(19) 



Since this is just retracing the path in which we obtained 



The result is shown by the dash-dot line in Fig. 4. Evi- 
dently, the contribution from the thermal-shower recom- 
bination is much more important than that from frag- 
mentation in the range of pt shown. Despite the fact 
that T(fcr) is lower than L(/ct) for fc^ > 1.5 GeV/c, its 
dominance at fcy < 1.5 GeV/c is enough to result in the 
TS recombination to dominate over the SS recombina- 
tion for all pt < 8 GeV/c. This example demonstrates 
the necessity of knowing the SPD's in a jet, since K(x\) is 
used in Eq. I|17|) . If SS recombination is the only impor- 
tant contribution as in pp collisions, then fragmentation 
as in Eq. (|19fl is all that is needed, and the search for 
SPD's plays no crucial role. In realistic problems the 
hard-parton momentum fc™ ax has to be integrated over 
the weight of the jet cross section. However, for our il- 
lustrative purpose here, that is beside the point. 

Our next example is the study of the dihadron distri- 
bution in a jet. We need only carry out the investiga- 
tion here for a jet in vacuum, since the replacement of a 
shower parton by a thermal parton for a jet in a medium 
is trivial, having seen how that is done in the replace- 
ment of Eq. (|18fl by (|17fl in the case of the single-particle 
distribution. Consider the joint distribution of two tt + in 
a jet initiated by a hard u quark. As we shall work in the 
momentum fraction variables, the value of the momen- 
tum of the initiating u quark is irrelevant, except that it 
should be high. Let Xi and denote the momentum 
fractions of the two 7r + , and Xi denotes that of the ith 
parton, i = 1, • • • ,4. Then, since only one u quark can 
be valence, the other three quarks being in the sea, we 
have one K, three L, and two R functions. Combina- 
torial complications arise when we impose the condition 
that J2i ij < 1 for i = 1, 2, 3, 4. There are two methods 
to keep the accounting of the different orderings of the 
four Xi. 



Method 1. 

Let one ordering be 



SPD(x\X2X^X4) = K(x\)L 



X-2 



1 — X\ 



L 



x 3 



1 — X\ 



x 4 



1 — xi — xi — x 3 



(20) 



There are 4! ways to rearrange the four Xi in all orders. However, they are to be convoluted with i? Tr (a;i, xi, X\), 
which is symmetric in x\ <-> X2, and similarly with R n (xs, X4, X2). Thus there are 4!/2!2! independent terms. Since 
K can appear at any one of the four positions in Eq. I|2U|I . we have altogether 24 terms. Thus we have 



X X X. 



dN n + n + 
dX x dX 2 



Wi 1 J L p 



R v (xi,x 2 ,X 1 )R 7r (x 3 , X4, Xi) 



(21) 
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where symbolizes the permutation of all Xi and summing over all four positions of K, but eliminating redundant 
terms that are symmetric under the interchanges of X\ x 2 and x% <-> X4. 
Method 2. 

Let us fix the ordering in Eq. (|2UI) but permute the contributing Xi to X\ and X 2 . There are six arrangements of Xi 
and Xj in R n (xi, Xj, Xi)R K (xi> , xy , X 2 ), while counting in xy and xy is unnecessary. Let us denote the summation 
over them by X)o- Thus we have 



X\X 2 



dX x dX 2 



n 



dxi 



- ^2 SPD(x 1 X 2 X 3 X4) 



K 



^y^i^ (x< , , Xi )i? 7r (av , a;j/ , X 2 ) 
Q 



(22) 



where denotes summing over the four positions of K. Equation (1221 is equivalent to (|21(1 . 

It should be noted that not all terms in these equations can be expressed in the form factorizable FF's. One example 
that can is 



n 



dxj 



xL 



X3 



1 — X\ — X2 



K{ Xl )L 



L 



■r-2 



1 — X\ 

Xi 



+ L( Xl )K 



x 2 



1 — X\ — X2 — X3 



1 — X\ 

R 7r {x 3 ,x 4 ,X 2 



Ri T (xi,x 2 ,X 1 ) 



DZ (X 1 )D- S {X 2 /{l-X 1 )). 



(23) 



Because of the presence of terms that cannot be written 
in factorizable form, the two-particle distribution cannot 
be adequately represented by the FF's only. 



T3 



1 r 


— X i 


=0.4 




- X i 


=0.5 
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0.1 0.2 v 0.3 0.4 0.5 



FIG. 5: Two 7r + correlated distribution in a u-quark initiated 
jet. 

Using the SPD's obtained in the previous section, we 
get the results shown in Fig. 5, which exhibits the X 2 dis- 
tribution for four fixed values of X\ . This type of corre- 
lation in parton fragmentation has never been calculated 
before. Although the shapes of the X 2 distributions look 
similar in the log scale in Fig. 5, there is significant at- 
tenuation as X 2 — ► 1 — X\ for each value of X\. Thus 
the effective slope becomes steeper for larger X\ . Recent 
experiments at RHIC have begun to measure the distri- 
bution of particles associated with triggers restricted to 



a small interval. The extension of our calculation here to 
such problems will need the input of jet cross sections for 
all hard partons in heavy-ion collisions and the partici- 
pation of thermal partons in the recombination. Here we 
only demonstrate the utility of the SPD's in the study of 
dihadron correlation. 



VI. CONCLUSION 

We have described the fragmentation process in the 
framework of recombination. The shower parton dis- 
tributions obtained are shown to be useful in problems 
where the knowledge of the fragmentation functions alone 
is not sufficient to provide answers to questions concern- 
ing the interaction between a jet and its surrounding 
medium or between particles within a jet. Such ques- 
tions arise mainly in nuclear collisions at high energies. 

In our view the basic hadronization process is recom- 
bination, even for fragmentation in vacuum. Since the 
recombination process can only be formulated in the 
framework of a model, the shower parton distributions 
obtained are indeed model dependent. That is a price 
that must be paid for the study of hadrons produced at 
intermediate pt where the interaction between soft and 
semi-hard partons cannot be ignored, and where per- 
turbative QCD is not reliable. Once recombination is 
adopted for treating hadronization in that pr range, the 
extension to higher pt can remain in the recombination 
framework, since the fragmentation process is recovered 
by the recombination of two shower partons. For hadron 
production in heavy-ion collisions at super high energies, 
such as at LHC, then the high density of hard partons 
produced will require the consideration of recombination 
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of hard partons from overlapping jets. Thus it is sensible 
to remain in the recombination mode for all pt- 

We have shown in this paper how the SPD's can be de- 
termined from the FF's. Although we have determined 
the SPD's at only one value of Q 2 for the FF's, it is clear 
that the same procedure can be followed for other value 
of Q 2 . The formal description of how the Q 2 dependences 
of the FF's can be transferred to the Q 2 dependences of 
the SPD's is a problem that is worth dedicated attention. 
While the numerical accuracy of the SPD's obtained here 
can still be improved, especially at lower Q 2 , for the pur- 
pose of phenomenological applications the availability of 
the parametrizations given in Table II is far more im- 
portant than not taking into account at all the shower 
partons and their interactions with the medium in the 
environment. The Q 2 evolution of the SPD's may have 



to undergo a long process of investigatory evolution of its 
own just as what has happened to the FF's. That can 
proceed in parallel to the rich phenomenology that can 
now be pursued in the application of the role of shower 
partons to heavy-ion collisions. 
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